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Abstract: 

The pseudo SU(3) model is used to study the double beta decay of ^^'^Nd to the ground 
and excited states of ^^^Sm. Low lying collective excitations of ^^^Sm and its BE(2) intensities 
are well reproduced. Expressions for the two neutrino double beta decay to excited states are 
developed and used to describe the decay of ^^^Nd. The existence of selection rules which 
strongly restrict the decay is discussed. 
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Introduction 



The neutrinoless double beta decay (PPou), undetected up to now, provides the more strin- 
gent limits to the Majorana mass of the neutrino < m,^^ >< l.ley||l|. Its detection would 
imply an indisputable evidence of physics beyond the standard model and would be useful in 
order to select Grand Unification Theories]^. 

Theoretical nuclear matrix elements are needed to convert experimental half-life limits, 
which are available for many /3/3-unstable isotopes0, into constrains for particle physics param- 
eters such as the effective Majorana mass of the neutrino and the contribution of right-handed 
currents to the weak interactions. Thus, these matrix elements are essential to understand the 
underlying physics. 

The two neutrino mode of the double beta decay {PP2u) is allowed as a second order process 
in the standard model. It has been detected in nine nuclei and has served as a test of a variety 
of nuclear models. The calculation of the and /9/?oi/ matrix elements requires different 
theoretical methods. Therefore a successful prediction of the former cannot be considered a 
rigorous test of the latter, but gives some confidence. However it is the best available proof we 
can impose to a nuclear model used to predict the PPqu matrix elements. 

Many experimental groups have reported measurements of (3(321, processes]^. Nearly for 
all the cases the ground state (g.s.) to ground state (0"*" 0^) decay was investigated. In 
direct-counting experiments the analysis of the sum-energy spectrum of the emitted electrons 
allowed the identification of the different /5/?-decay modes[|]. 

Recently the possibility of detecting PP decay into excited states of the daughter nucleus 
by measuring the gamma radiation has exerted some attraction among experimentalists. This 
is due to the fact that phase space integrals scale as the energy available for the decay and 
consequently decrease for excited states. In the case of the decay to a first excited 2"*" state the 
phase space factor contains terms which are antisymmetric in the energies of the two outgoing 
electrons and antineutrinos, resulting in a large reduction of the corresponding integral. It 
makes such transitions very difficult to observe]^, 0. Also they are inclusive experiments which 
cannot distinguish between the different (3(3 decay modes. On the other hand the detection of 
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a gamma ray gives a much clearer signal than a continuous electron spectrum as in the case of 
the g.s. — > g.s. decay. 

The pioneer work of Bellotti et. al^ determined lower limits to the I3j3 decays to excited 
states of six nuclei. Lower limits for the decay of ^"^^Te and ^^^Te to the first 2"^ state of the 
corresponding Xenon isotopes have been reported^]. The (3(3 decay of ^^Ge to excited states 
of "^^Se was studied looking for the detection of one or two photons in coincidence with two 
electrons|P, |lOl- The half-lives for the (i(i2u decay from the g.s. of ^°°Mo, ^^Zr and ^^^Nd to 



the first excited 0"^ state of the daughter nuclei have been estimated for the first time in |11 



assuming that the nuclear matrix elements of the transitions to the g.s. and excited 0''" state are 
the same. Experimental studies of the decay to excited states looking for the 7 ray signature 
have been performed for ^^Zr[||], i^^Cdl, |13|, H), and i°°Mo[|13], [l^. The detection of 



the (3(32u to the first excited O"*" state was reported for the first time in |[T^, The factibility 
of studying the (3(32v to excited states in ^^^Nd has been discussed in recent yearsQ, [11], |1^, |19|. 
In [|12[ preliminary results were reported. 

Theoretical analysis of the /5/?2i/ to excited states have been performed in the context of 
the QRPA formalism for ^^^Mo^, ^^^Xe^, |2| and ™Ge0, gSf, "^Crfg and also for 
^"^Se,^^^ Pd and ^^^'-"^^^Tel^BI. When the first excited O'*' state was studied it was assumed that 
it is a member of a two quadrupole-phonon triplet. The QRPA calculation for ^''^Mo exhibits 
an overestimation of the amplitude of the (3(32u decay to this excited state when the decay to 
the ground state is reproduced[^, For the case of ^^^CrfQ the matrix element of this decay 
is five times greater than that associated with the decay to the ground state. 

Although there is no reported calculation of the (3(32v of ^^^Nd to excited states of ^^^Sm, 



this nucleus was mentioned as a suitable candidate for this decay. In this conclusion has 
been reached assuming that the nuclear matrix elements for the fifi2u to the ground state and 
the first excited 0"^ state are equal. In Q it was speculated that if the matrix element for the 
/?/?2v of ^^^Nd would show a similar enhancement over that of the g.s. decay as found for ^^^Cd, 
the decay rate into this excited level could even exceed that of the g.s. decay. 

In the present paper we perform an analysis of the (3(32u decay of ^^^Nd to excited 0^ and 
2+ states of ^^^Sm using the pseudo SU(3) formalism. This theoretical model is well suited to 
describe the collective spectra. Under this scheme the first 2"^ and 4"^ states are members of 



3 



the g.s. rotational band, while the second excited 0^ and 2+ states are members of an excited 
rotational band. The third excited 0^ state is the head of another rotational band. We will 
show that whithin the pseudo SU(3) model the (3(32u decay to the first excited O'^ state will be 
cancelled while the other decays will be strongly suppressed. The necessary formalism to study 
the PP2u decay to 2+ states is also developed. 

In Section 2 the pseudo SU(3) formalism and the model Hamiltonian are briefly reviewed. In 
Section 3 the summation method is used to obtain the (3(32u matrix elements for the decay to the 
2"*" states. Section 4 contains the explicit formulae needed to evaluate the l3l32u matrix elements 
in the pseudo SU(3) scheme. The nuclear structure analysis of ^^^Sm is given in Section 5. In 
Section 6 the l3l32u nuclear matrix elements and half-lives are presented. Conclusions are drawn 
in the last Section. 



The pseudo SU(3) formalism 

In order to obtain a microscopical description of the low lying energy states of ^^^Nd and 
^^^Sm we will use the pseudo SU(3) model which successfully describes collective excitations in 
rare earth nuclei and actinides|2^ as well as the g.s. ^ g.s. (3^2^ and PjSoiy decays of six heavy 
deformed nuclei^ ^ §D|. 

In the pseudo SU{3) shell model coupling scheme PT| normal parity orbitals {ri,l,j) are 



identified with orbitals of a harmonic oscillator of one quanta less fj = rj—l. This set of orbitals 
with j = j = I + s, pseudo spin 5 = 1/2 and pseudo orbital angular momentum / define the so 
called pseudo space. Recently it was found an analytic expression for the transformation of the 
normal parity orbitals to the pseudo space |32|. Applying this transformation to the spherical 



Nilsson Hamiltonian it can be shown explicitly that the strength of the pseudo spin orbit 
interaction is almost zero for heavy nuclei and the orbitals j = I ±1/2 are nearly degenerate. For 
configurations of identical particles occupying a single j orbital of abnormal parity a convenient 
characterization of states is made by means of the seniority coupling scheme. 

The many particle states of nucleons in a given shell 7]^, a = u or tt, can be defined by 
the totally antisymmetric irreducible representations } and of unitary groups. The 
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dimensions of the normal (A^) parity space is = {fja + l)(^7a + 2) and that of the unique (A) 
space is = 21]^ + 4:, with the constraint Ua = n^ + n^. Proton and neutron states are coupled 
to angular momentum and in both the normal and unique parity sectors, respectively. 
The wave function of the many-particle state with angular momentum J and projection M is 
expressed as a direct product of the normal and unique parity ones, as: 



\JM >-- 



E [|J^>®|J^ >]i, 



We are interested in the low energy states of ^^^Nd and ^^^Sm which have J = 0, 2, 3, 4, 6. 

For even-even heavy nuclei it has been shown that if the residual neutron-proton interaction 
is of the quadrupole type, regardless of the interaction in the proton and neutron spaces, the 
most important normal parity configurations are those with highest spatial symmetry {fa} = 
|2"a/2j|g^ This statement is valid for yrast states below the backbending region. This implies 
that = Sj, = 0, i.e. only pseudo spin zero configurations are considered. 

Additionally in the abnormal parity space only seniority zero configurations are taken into 
account. This simplification implies that = = 0. This is a very strong assumption 
quite useful in order to simplify the calculations. Its effects upon the present calculation are 
discussed below. 

The double beta decay, when described in the pseudo SU(3) scheme, is strongly dependent 
on the occupation numbers for protons and neutrons in the normal and abnormal parity states 

These numbers are determined filling the Nilsson levels from below. 



N N A A 



28, 29 



as discussed in p8|, |29[. In particular the j3j32u decay is allowed only if it fulfils the following 
relationships 



N 



+ 2 



n. 



N 



N 



<.-2 



(2) 



For ^^^Nd it was assumed a deformation P ^ 0.28 |Q. We have obtained the occupation 
numbers 



n 



N 



n 



N 



(3) 
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For ^^^Sm it is found that (3 ^ 0.19|Q but we were forced to select a higher deformation to 



satisfy relations ( According to this higher deformation is more appropriate for ^^'^Sm 
and is related with some departure from a rotational behavior in the ground state band of 
^^°S'm. The selected occupation numbers for ^^^Sm are 

n;f = 6 , = 6 , n;^ = 2 , = 4 . (4) 

In order to analyze the spectra and transitions amplitudes of ^^°S'm we have selected the 
standard version of the pseudo SU(3) Hamiltonian [^] . It is constructed by a spherical central 
potential, a quadrupole-quadrupole interaction and a residual force. The latter allows the fine 
tuning of low lying spectral features like K band splitting and the effective moments of inertia. 
The hamiltonian looks 

H = Y.Ha- IxQ^ ■ Q"Ci + C2 L'. (5) 

The spherical Nilsson hamiltonian which describe the single-particle motion of neutrons 
[a = v) or protons (a = vr) is: 

Ha = ^hu} irias + I - 2A;a L„s ■ S^s - Kl^aL'iX - = ^ CsaCtL^sa (6) 
s ^ ^ s,a 

where rj = fj + 1 denotes the harmonic oscillator number operator and hu determines the size 
of the shell. A constant term (Vyr) is included which represents the depth of the neutron 
(proton) potential well. In (6) the second quantization representation of Ha is given, ega being 
the single-particle energies. 

The quadrupole operator = {l-ns + Ivs) ^'^^^ o^^Y within a shell and do not mix 
different shells. The residual interaction, K^, is a linear combination of L^, and X4, defined 
as 

Xs= T.yUQ'i^L, ^ (7) 
= Y.i,j,k LiQijQ'jkLk 

They are rotational invariant and scalar operators built by generators of the algebra of SU{3) ^7, 



9[ and and Q'^j are cartesian forms of the total angular momentum and the quadrupole 
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operators, respectively. The K is interpreted to be the third component of the total angular 
momentum on an intrinsic body-fixed symmetry axis of the system, which is given by 



= (A1A2L2 + A3X3 + X,)/{2Xl + A1A2) (8) 

with the parameters Aj denoting the eigenvalues of the mass quadrupole operator, which are 
related to the SU{3) labels (A, /i) through the expressions 

Ai = i(/i-A), A2 = -i(A + 2/x + 3), A3 = ^(2A + ^ + 3). (9) 

Although the quantum number K used to define the orthonormalised basis is not the same 
as the Elliot k, the states studied in the present work satisfy quite accurately the relationships: 

K'^\K = 1>=0, K'^\K = 2>^A\K = 2> . (10) 

It would be possible to add the the hamiltonian terms which distinguish different irreps. 
For the sake of simplicity we keep this simplest version. 

With the occupation numbers determined in Eq and and the hamiltonian the 



wave function of the deformed ground state of °"Nd can be written 28, 29 



|i50^rf,0+) = |0+) = I = = 0; (^13/2)^, = = >a 

I {16}423},(12,0V; {16}423},(18,0),; 1(30, 0)ir = IJ = M = >;v , 

(11) 

and the deformed low energy states of ^^^Sm are described by the wave functions 



rSm, J+) = 14+) = I ihn/2t, = = 0; (213/2)^, = = >a 

I {16},{23},(12,0)^; {14},{22},(12,2),; 1(A, J M >jv , ^ ' 

where J+ denotes a state with angular momentum J, positive parity and associated with the 
SU(3) irrep (A, /x)^-. In this approach we are assuming that the first 0+, 2+, 4+ states of ^^^Sm are 
the low energy sector of a rotational band described by the normal (A, ^j)g.s. = (24, 2) strong 
coupled pseudo SU(3) irrep, the second 0"'",2+ states belong to a second rotational band with 
(A,yu)i = (20,4), and the third 0+ state is the head of another rotational band described by 
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the pseudo SU(3) irrep (A,/i)2 = (22,0). We will discuss also a gamma band associated with 
(A, fJ')g.s. with K=2. 



The PP2U decay to excited 2+ states 



The inverse half life of the two neutrino mode of the /3/3-decay can be cast in the form[^ 



1/2/ 
'2u 



-1 



r - j:) = I MM) \' . (13) 



where G2v{J^) are kinematical factors. They depend on Ej.,j = ^[Q/313 — E{J^)] + rrieC^ which 
is the half of total energy released. The nuclear matrix element is 

MM) - MfJ( J+) = E " " '"M'f " " (14) 

with the Gamow- Teller operator F expressed as 

^m = J2 ^mst; =J2('i^^'')H^l^^-j^^%J.i,jfm m = 1, 0, -1. (15) 

S TT U 

The energy denominator is fi^ = Ej^^ + Ej^ — Ei and it contains the intermediate E^ and 
initial Ei energies. The kets |1^) denote intermediate states. 

The mathematical expressions needed to evaluate the nuclear matrix elements of the g.s. —>■ 



g.s. f3f32u decay in the pseudo SU(3) model were developed recently[^, The same formulae 
describe the decay to the first excited 0+ state by replacing the values of the strong coupled 
irrep (A, /u)g.s. of Eq. (|12|) with those corresponding to excited bands. 

We will concentrate first on the derivation of the matrix element M2,y(2+) to the 2"^ excited 
states. The formulae for this decay resembles that of the decay to the O"*" states but the energy 
denominator is up to the third power. Being in general this energy of the order of 10 MeV 
this power implies a factor 100 of suppression for this matrix element]^, The previous 

equation is rearranged as 

Mf,^(2+) = v^5](l/il/i'|2m) E /i^^(2>| T,, \ l%m,) (l+m^ | \ 0+) . (16) 

li/i' Nm\ 
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Using 



1 

2 dEl^ 



;i7) 



and the summation method described inp9|, ^ it is possible to rewrite the second sum as 



S |J^W(2>|r^'| |rj 0+)} = 

ial^(2>| g ^r,,[i/, [H,..., [H, T,] . . .](A-t-e.)|o+^ 



The two body terms of the Hamiltonian (^) commutes with the Gamow- Teller operator 
(P^), thus the above multiple commutators are easy to evaluate. We obtain [P^ 



[H,...[H,T,]]... 



Enu z/)[at ® a^]i^{e^ - ej^ 



(19) 



where vr = (//tt, ^tt, j-Tr) and = (rj^Juyju)- Returning with this expression to the original 
formula, Eq. (13), resumming the infinite series and recoupling the Gamow- Teller operators it 
is found 



1 a2 

2 dE'i 



[4/ a^'Y 



2m 



Of) 



2m 



(20) 

As it was shown in |^ the expression for the nuclear matrix element of the (3(32u decay to 
a 0"*" state is similar to (PDj), with a different power in the denominator. Equations ( [2D| ) and 
(4.10) of can be expressed in a compact form as: 



Mf7(J+ 



-(J+m| 



J?™ 



0+) 



1 



-(J+m|T^'"(7rz/,7rV')| 0+) 



(21) 



For practical purposes the tensor T"^™(7rz/, vrV) was implicitly defined in the above equation. 
The 73 for the J = case comes from the relation between the scalar product of two vectors 
and their coupling to angular momentum zero. 
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The matrix elements 



We want to evaluate the nuclear matrix element ( pi]) for the l3P2u decay of the ground 
state of ^^^Nd, Eq. ([llD , to the ground and excited states of ^^^Sm which are described by 



the wave functions of Eq. ([T2|) . Each Gamow- Teller operator ( p!5|) annihilates a proton and 
creates a neutron in the same oscillator shell and with the same orbital angular momentum. 
In the case of the (3(32u of ^^^Nd it means that the operator annihilates two neutrons in the 
pseudo shell rji, = 5 and creates two protons in the abnormal orbit hii/2- As a consequence 
the only orbitals which in the model space can be connected through the PP2U decay are those 
satisfying 77^ = 771^ = 77, that implies K = lu = Jv = V " k jn = V + \- These are the 
selection rules described by relations (2) concerning the change in occupation numbers. Under 
this restrictions only one term in the sum over configurations vrz/, vrV survives and thus the 
nuclear matrix element M21, (|2T|) can be written as 



= -^,{J^\T'-\nu,nu)\Ot), (22) 

where the energy denominator is determined demanding that the Isobaric Analog State in the 
intermediate odd-odd nucleus is an eigenstate of the Hamiltonian (|^). Its excitation energy is 
equal to the difference in Coulomb energies A^. Their expressions are[2£] 

= Ej„ + e(?7^, = jy + 1) - e(?7j., l^^j^) = Ej„ - hujk^2j^ + Aq. 

(23) 

Ac = ^[2Z+l- 0.76((Z + 1)4/3 _ z^/^)]MeV. 



As it was discussed in |^ in the context of the g.s. g.s. (3(32u decay Eq. ( ^2]) has no 
free parameters, being the denominator (|23| ) a well defined quantity. The reduction to only one 
term comes as a consequence of the restricted Hilbert proton and neutron spaces of the model. 
The initial and final ground states are strongly correlated with a very rich structure in terms 
of their shell model components. 

Since the Hilbert space has been divided in their normal and unique parity components we 
need to rearrange the creation an annihilation operators in the same way, i.e. 
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\ai (g) a,y 



[al (g) a,yY 



JM 



[40 4]-^- ® [h^®hyY" 



JM 



(24) 



The is the unitary (Jahn-Hope) 9-j recouphng coefficient!^. Introducing this ex- 

pression in (p^) together with the exphcit form of the wave functions ([l^) and (p!2D we obtain 



ivr ju 1 



JivJiy \ J IT 

^ TV ^ U 

< {hu/2)l, = = 0|[4 ® alY^l {hu/2)t, = M^ = 0>® 
< (12,0)^; (12, 2),; 1{X, fi)^K = lJM\[d, (g) d^Y^ 

I (12, 0)^; (18, 0)^; 1(30, 0)K =1J = M = >y^^ 

(25) 

The matrix element (^) vanishes unless a pair of protons coupled to total angular momen- 
tum zero is created and two neutrons of the normal parity space coupled to pseudo orbital 
angular momentum L = J and pseudo spin equal to zero are annihilated. The above sum is 
thus restricted to = 0, Ju = J- 

The operators in the normal space must be recoupled from the jj— to the LS*— coupling 
scheme. The result is 



[ 1 



Ju 



iLS 

(0V)1^\JM 



(26) 



LS 



L S J 



The low energy levels are assumed to have pseudo spin S = 0, a fact that again simplifies 
the evaluation of the above sum by imposing L = J. 

Using that = + IJtt = lu = V the reduced matrix elements cr(7rz/) read 



l7]{r] + 1) 



(27) 



3(2r7 + l)' 

In the seniority zero approximation the two particle transfer matrix element of the unique 



sector of Eq. (22) is evaluated by using the quasispin formalism and gives p8|, p9 | 
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^ TV 'TV TV I L TT ^ TTJ I TT ' TV TV 



>= 



(n^ + 2)(r/ + l-n^/2) 
?7 + 1 



1/2 



(2^ 



The evaluation of the matrix elements in the normal space of Eq. p5|) is performed by 
using SU (3) Racah calculus to decouple the proton and neutron normal irreps, and expanding 
the annihilation operators of Eq. (pUf ) in their SU(3) tensorial components. The final result is: 



MfJ(J+) = a(J) bin^) £jy^'^ 

E < mil, {Ofi)ll\\{Xofio)KoJ >i E < (30,0)1 0, {Xofio)KoJ\\{Xfi)AJ >^ 

i\oiio)Ko P 



E 

p' 



(12,0) (0,0) (12,0) 1 

(18,0) (Ao/io) (12,2) p' 

(30,0) (Ao/io) (A/i). p 
111 



(29) 



<(12,2)|||[ao^),iao^),i](^°^-Mll(18,0) >p. 



In the above formula < 



> denotes the SU(3) Clebsch-Gordan coefficients!^, the 



symbol [...] represents a 9 — A/i recoupling coefficient |]38[| , and < .. 



.. > is the triple 



reduced matrix elements [^. The energy denominator was defined in (23), and 



a(0) 



a(2) 



5>7(r?-l))2r;-3) 
3(27?+l) 



1/2 



(30) 



The rotational spectrum of ^^^Sm 

The four parameters of the pseudo SU(3) hamiltonian (|) were fitted to reproduce the first 
2+ states in ^^"S'mas it was done in Ref.fl^. Their values are 

X = 3A7eV , Ci = 2l5keV , (2 = 50AkeV . (31) 

The right hand side of Fig. 1 exhibits nine of the lowest energy states which have been 
observed in ^^^Sm, grouped in rotational bands. Angular momentum and parity are given 
for each level. The left hand side of Fig. 1 shows the calculated spectrum together with the 
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associated irreps. The gamma band is identified with K=2. The general trend is well reproduced 
but the experimental g.s. band does not show the rotational structure which is exhibited by the 
calculated one. This departure from an exact rotational behaviour was mentioned in Section 2 
and it can be associated with the relatively small deformation reported for ^^^Sm. In this mass 
region the deformation suddenly jumps for ^^'^Sm to a value very similar to the deformation of 
^^^Nd. The gamma band is present in the g.s. irrep because both A and /i are different from 
zero and is well fitted. 

The excited O"*" states are the head of other rotational bands. The predicted energy gap 
between them is 125 eV while the experimental one is 454 eV. These numbers suggest that we 
have not a clear identification of these excited states. This fact has relevance in the study of 
the /3/32ydecay to these states. 



The BE(2) transition intensities were evaluated using the effective quadrupole operator p7 



Qo = el"Q. + et"Q, , el" = e + epoi , et" = e^oi (32) 

with Cpoi = 0.93e. The seniority zero condition imposed to the nucleons in abnormal parity 
orbitals inhibits them to participate in collective excitations. This restriction forces a slightly 
large value for the polarisation charge. A similar effect was found in a BE(2) study of rare 
earth and actinide nuclei p7[] . 

The BE(2) intensities of the transition from the first and second 2"*" to the ground state, and 
from the first to the state are shown in Table 1 and compared with their experimental 
values, in Weiskoff units (W.u.). The agreement is good except for the case of the transition 
from the second excited 2"*" state to the ground state which fails in a factor four. 



The /?/?2.. decay of ^^^Nd 



In this section we study the two neutrino mode of the double beta decay (3(32u of -^^^Nd into 
the ground state, the first excited 2"*" and the first and second excited O"*" states of ^^^Sm. 

In Table 2 the matrix elements, energy denominators, phase space integrals and predicted 
half-lives for the /3/32iy decay of ^^^Nd to the ground state, the first 2"*" and the first and second 
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excited 0"*^ states of ^^^Sm are presented. The matrix elements (26) are given in units of 
)-(J+i)_ Phase space integrals for the decay to 0"*" states were evaluated following the 
prescriptions given in|^ with qa/ Qv = 1-0 and the kinematical factor for the decay to the 2"*" 
state was taken from |^ renormalized by the above mentioned value of the axial vector coupling 
constant. It must be mentioned that these phase space factors differ in about 10% with those 



used in [psj P9|] where a different renormalization procedure was used. 



As it was mentioned in |2^ the predicted half live for the (3f32v to the ground state of 



^^'^Sm is in reasonable agreement with the experimental data, which vary between 9 and 17x 10^^ 
years!, 

The (3(32v decay to the first excited 0^ state is forbidden. In this model this is imposed by 
the fulfilment of an exact selection rule. It can be understood by realizing that the pair of 
annihilation operators 5(04)1) when expanded in their SU(3) components, have the couplings 
(0,4) ® (0,4) = (0,8), (2, 4) containing a L = state. But acting over the ^^^Nd g.s. irrep 
(30,0) they cannot couple to the irrep (20,4) which we associated with the first excited O"*" state. 
In other words the transition between members of these particular irreps are forbidden. 

The decay to the second excited 0"^ state is allowed but strongly cancelled. The reduction 
of the matrix element by a factor ten, in comparison with that associated with the decay to the 
g.s., is partly related with the fact that though the coupling (30, 0) ® (0, 8) = (22, 0) is allowed 
the coupling with (2,4) is forbidden. The predicted half-life is four orders of magnitude larger 
than that of the decay to the g.s. 

The (3l32v decay to the 2"^ state is inhibited by the /i^ dependence of the matrix element as 
it is discussed in Section 3. The matrix element of the /3/32;ydecay to the first excited 2+ state 
M2v{2~^ g ) is three orders of magnitude lesser than the matrix element of the decay to the g.s. 
A^2.(0+J. 

The present results contradicts those previously published |Q, ^ were it was speculated that 
the j3j32u decay of ^^^Nd to the first excited O"*" state of ^^^Sm could have a similar intensity 
of that to the g.s. We find that in the present formalism this decay is forbidden. If we select 
different occupation numbers for both ^^^Nd and ^^^Sm, taken the deformation of the latter 
nucleus instead of that of the former, we find very similar results for the decay to the g.s. 
and the first 2^ state, but the matrix elements of the decay to the first and second 0"*" states 
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becomes interchanged with essentially the same values. Considering the difference in the phase 
space integrals we predict a half live of the order 10^^ years for the decay to the first excited 
0"*" state and the decay to the second excited one becomes forbidden. 

The above discussed reduction of the matrix element of the decay to the excited 0"*" 
state as compared with the decay to the g.s. is not a general result of the pseudo SU(3) 
scheme. A recent analysis of the case of ^°°Mo^^ shows that both matrix elements are very 
similar and that they are in agreement with the experimental information. In conclusion, 
the appearance of selection rules which can produce the suppression of the matrix elements 
governing a (3(32u transition is a consequence of the details of the irreps involved. 

The pseudo SU(3) model uses a quite restrictive Hilbert space. The model could be improved 
by incorporating mixing between different irreps, via pairing by example ||44]| . Also other active 
shells can be taken into account in the symplectic extensionj^^. In both cases the selection 
rules that impose such strong restrictions to the decays of some nuclei can be superseded. 
However if the main part of the wave function is well represented by the pseudo SU(3) model 
those forbidden decays will have, in the better case, matrix elements that will be no greater 
than 20% of the allowed ones, resulting in at least one order of magnitude cancellation in 
the half-life. In any case these results should be taken into account in the design of future 
experiments. 



Conclusions 



In the present paper we have studied the (3(32u decay mode of ^^^Nd to the ground and 
excited states of ^^^Sm. The transitions have been analyzed in the context of the pseudo SU(3) 
model. The experimental spectrum of the g.s. rotational band of ^^^Sm was reproduced as well 
as the measured half-life of the l3(32u decay to the g.s. but the excited rotational bands were 
not so well reproduced. The PP2U decay to the first excited 0"^ state was found forbidden in 
the model and the decay to the second excited 0^ state has a half-life four orders of magnitude 
greater that that to the g.s.. The decay to the 2+ state is strongly inhibited due to the energy 
dependence of the matrix elements M2jy(2+), two powers greater than that of the matrix element 
M2,(0+). 
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It is expected that improving the model would remove the exact selection rules which forbid 
some decays. In any case, if the pseudo SU(3) wave functions are a good representation of the 
low-lying energy states of ^^^Nd and ^^^Sm, they will remain inhibited. 
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Table captions 



Table 1 The BE(2) intensities, in Weiskoff units (W.u.), of the transition from the first and 
second 2+ to the ground state, and from the first 4+ to the 2+ state of ^^^Sm are shown. The 
calculated and experimental values are exhibited in the second and third column, respectively. 

Table 2 The dimensionless /3/32i/ matrix elements, energy denominators, phase space inte- 
grals and predicted half-hves for the decay of ^^^Nd to the ground state, the first 2+ and the 
first and second excited 0+ states of ^^^Sm. 



Figure captions 



Figure 1 Spectrum of the low-lying states of ^^^Sm. The levels are grouped in rotational 
bands and they are labeled by angular momentum and parity. The right hand side contains 
the experimentally determined levels. In the left hand side the calculated spectrum is exhibited 
with the associated irreps at the bottom. 
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Table 1 



transition theory experiment 



2t- 




55.7 


55.8 


2t- 




0.48 


2.0 


4^- 


-2+ 


78.5 


112.0 



Table 2 







Sj^iMeV] 


G,^{J^)[yr-'] 




0+- 




.0549 


12.20 


4.94 X 10-^^ 


6.73 X 10^^ 


0+- 


-o+(i) 





11.83 


5.83 X 10"^® 


oo 


0+- 


- 0+(2) 


.00499 


11.58 


9.33 X 10-19 


4.31 X 10^2 


0+- 


-^2+ 


5.38 X 10-5 


12.04 


4.78 X 10-1^ 


7.21 X 10^4 
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